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Evidence for filamentary superconductivity nucleated at antiphase domain walls in 
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Resistivity, magnetization and microscopic '^^As nuclear magnetic resonance (NMR) measure- 
ments in the antiferromagnetically ordered state of the iron-based superconductor parent material 
CaFe2As2 exhibit anomalous features that are consistent with the collective freezing of domain walls. 
Below T* ~ 10 K, the resistivity exhibits a peak and downturn, the bulk magnetization exhibits a 
sharp increase, and '^^As NMR measurements reveal the presence of slow fluctuations of the hyper- 
fine field. These features in both the charge and spin response are strongly field dependent, are fully 
suppressed by H* ~ 15 T, and suggest the presence of filamentary superconductivity nucleated at 
the antiphase domain walls in this material. 



A. Introduction 

The interplay among competing ground states of cor- 
related electron systems can give rise to a rich spectrum 
of emergent behavior. The iron-based superconductors 
are particularly noteworthy, and have attracted extensive 
interest since the discovery of La[Oi_a;Fa;]FeAs in 2008.^ 
Like the high temperature superconducting cuprates, su- 
perconductivity (SC) in the iron arsenides emerges from 
an antiferromagnetic (AFM) parent state upon doping 
with excess charge carriers, and the superconducting 
pairing mechanism may be related to the AFM insta- 
bility of the parent state. ^"^ In both cases the electronic 
degrees of freedom condense into an unusual coexistence 
of both AFM and SC order parameters for intermedi- 
ate dopings.^'^ The nature of this coexistence is poorly 
understood. In general, a subdominant order parameter 
can emerge locally in regions where the dominant order 
vanishes or is suppressed.® Recent experiments in the iron 
arsenides suggest that SC and AFM order parameters are 
indeed spatially modulated on a microscopic scale. ^'^° 

Among the iron based superconductors, the AFe2As2 
(A = Ca, Sr and Ba) materials arc of particular interest 
because large single crystals of these oxygen-free com- 
pounds can be easily synthesized.^^ These materials un- 
dergo a tetragonal to orthorhombic transition followed 
by an AFM state upon cooling. ^^ Both chemical doping 
and applied pressure suppress the magnctostructural or- 
der of the parent compounds and give rise to SC, and 
the phase diagrams arc similar to those of other uncon- 
ventional superconductors. ^d3 q^j^^ CaFe2As2 system is 
noteworthy because SC is induced at only 0.4 GPa (non- 
hydrostatic), whereas the Sr and Ba materials require 2.8 
and 2.5 GPa, respectively. ^^"^^ 

In this paper we present evidence for coexisting fil- 
amentary SC and AFM in the undoped parent com- 
pound CaFe2As2, in which the SC order remains local- 
ized within AFM domain walls (DWs). A similar phe- 



nomenon has been observed in heavy fermion materials, ^^ 
and enhanced superfluid density has been observed at 
twin boundaries in Ba(Fei_2;Coa;)2As2."'^° CaFe2As2 is 
noteworthy, however, because the SC never achieves bulk 
long-range order, and surprisingly the filamentary SC ap- 
pears to be related to the presence of low frequency spin 
fluctuations. Measurements of the magnetotransport, 
magnetization, and nuclear magnetic resonance (NMR) 
reveal anomalies at T « 10 K that are suppressed with 
magnetic field in a manner consistent with the suppres- 
sion of bulk SC in doped samples. Similar anomalies 
in the AFM state were reported recently to correlate 
with the superconducting volume fraction in SrFe2As2.^® 
These features arc subtle in the CaFe2As2, but are mani- 
fest in both spin and charge probes. NMR spectra, spin- 
lattice (T]~ ) and spin-spin (^T^ ) relaxation measure- 
ments indicate slow dynamics and motional narrowing 
consistent with the freezing of mobile DWs.^^'^° The re- 
sistivity and magnetization both change abruptly at the 
same temperatures and fields and analysis of the mag- 
netotransport indicates the presence of weakly pinned 
superconducting filaments. This evidence suggests that 
the nucleation of finite superfluid density at DWs in 
CaFe2As2 drives the freezing of mobile antiphase do- 
mains. We thus find that this nominally pure stochio- 
metric material can spontaneously become electronically 
inhomogcneous at the nanoscale as a result of coupled lat- 
tice strain, antifcrromagnctism and superconductivity.^^ 



B. Experimental Details 

High quality single crystals of CaFe2As2 with typical 
dimension '^ 2x2x0.1 mm'^ were grown in Sn fiux by 
standard methods. ^^ Microprobc analysis indicates a Sn 
concentration < 600 ppm in the bulk. The in-plane resis- 
tivity Pah was measured using the electrical contact con- 
figuration of the flux transformer geometry. ^^ Six elec- 
trodes were fabricated on each sample by bonding Au 
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FIG. 1. (Color online) (a) T~^ vs. T for Ho 
♦ ,satellite),3T(^),6T(^)and9T(H) {H \\ c). (b) Tj ' vs 
Ho atSK {H \\c). (c) Tj"^ for H ± c at 9 T (solid green line 
is the like-spin second moment contribution.'^'^) (d) Fraction 
/ = 25/ A of signal intensity from DWs (see text for details) 
as a function of T. (e) Schematic of Fe spin orientations in 
plane indicating different types of domain walls: dotted blue 
lines are twin boundaries, dashed red are antiphase bound- 
aries along the a-axis, and dash-dot green lines are antiphase 
boundaries along the 6-axis. (f ) Spectra of the upper central 
transition (Hhyp \\ H || c) at several temperatures; the red 
lines are fits as described in the text. 



wires to the crystal with H20E epoxy paste and the cur- 
rent / was applied in the a6— plane. The NMR spectra 
and relaxation measurements were acquired on a single 
crystal at the upper of the two magnetically split central 
(/ = +i^-!>— i) resonances of the ^^As (/ = 3/2) using 
standard pulse sequences. 



C. Results and Discussion 

T^^ measurements in the AFM state of CaFe2As2 re- 
veal a small peak at 10 K that was attributed recently to 
the presence of slow spin fluctuations possibly associated 
with DW motion. ^^ To test this hypothesis, we measured 
Tj~ as a function of temperature in different applied field 
(Fig. 1(a), H II c) Surprisingly, we find that the intensity 
of this peak is suppressed by magnetic fields on the order 
of 10 T (Fig. 1(b), H II c). We also measured the spin- 
spin decoherence rate T2~^ versus T ior H ± c (Fig. 1(c)). 
In this configuration, T^^ = AujI^ + {l^ hl^^T,{T)f , 

where Aw^^p = 1.08 x 10^ sec~^ is the temperature inde- 
pendent second moment of the like-spin dipolar couplings 
among the ^^As for the central transition with H||[100] 
(solid line in Fig. 1(c)). ^■^'^'' Tg"^ clearly increases sharply 



below 25 K, revealing the presence of slow fluctuations 
of the hyperfine field in the ab plane. 

In the AFM ordered state, there arc two types of do- 
main boundaries that can give rise to slowly fluctuating 
hyperfine fields: twin boundaries and antiphase bound- 
aries (see Fig. 1(e)). ^° Twin boundaries emerge to relieve 
lattice strain and are immobile. ^^ Antiphase DWs have 
been detected well below Tn in BaFe2As2, and are likely 
to be mobile depending on the temperature.^^ Both types 
of DWs give rise to perturbations of the local hyperfine 
field and could be responsible for the slow dynamics we 
observe. ^^ The DWs give rise to a low frequency tail in 
the spectra that is clearly evident in Fig. 1(f). By using 
the known hyperfine couplings to the As, we model the 
resonance frequency as a function of position upon cross- 
ing such a DW as w(x) = ^Hq + wq tanh(a::/i5) where x 
is the position perpendicular to the DW, 5 is the DW 
width, ujq is the frequency arising from the hyperfine 
field, wo = 7i7hyp = 19.2 MHz, and -fHo = 65.4 MHz.^s 
We fit the spectra and extract the fraction / = 26/ A 
of the signal intensity arising from DWs, where A is the 
domain width. This portion is roughly independent of 
temperature down to 10 K, and then it decreases by ap- 
proximately 5% (Fig. 1(d)) and appears to be correlated 
with the slow dynamics observed in Ti and T2. Twin 
boundaries are most likely the source of the temperature 
independent contribution. Motional narrowing from mo- 
bile antiphase domain walls may be responsible for the 
the slow dynamics below 10 K, which may serve to sup- 
press 6 and hence /. In this case the correlation time 
'^c{T) > (jhi^yp)^^ ~ 10^^ s, where 7 is the gyromag- 
netic ratio and /ihyp ^ 1-4 kOe is the fluctuating com- 
ponent of the hyperfine field in the ab plane at the DW, 
which agrees with the calculated value at an antiphase 
boundary. 

We now turn to the bulk transport and magnetiza- 
tion experiments. Figure 2(a) depicts the temperature 
dependence of the resistivity pab a.t H = and 14 T 
{H II c), which exhibits a discontinuity at Tn =169 
K consistent with the reported phase transition. ^^ The 
magnetoresistance Ap/p = [p(14r) — p(0)]/p(0) (Fig. 
2(b)) is almost zero for T > T/v, but increases below 
T/v and shows a sharp upturn at 10 K. Figure 2(c) 
shows the zero-field-cooled (ZFC) and field-cooled (FC) 
magnetization M. The hysteresis below Tn suggests 
the presence of magnetic domains, and the low tem- 
perature increase in M may be due to free moments 
present in the DWs; we estimate the percentage of free 
moments to be ~7-18.7% for ordered moments of iron, 
/ie//=0.3-0.8 pb- Figure 2(d) depicts the T dependence 
of AM/M = [M(FC) - Af (ZFC)]/A/(ZFC). Below 10 
K the slope increases dramatically, indicating a collec- 
tive freezing of the magnetic domains. It is clear from 
these data that the charge transport is strongly coupled 
to the low temperature anomaly present in the NMR and 
magnetization data. 

In order to highlight the low temperature magneto- 
transport anomaly. Fig. 3(a) shows the low tempera- 
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FIG. 2. (Color online) (a) Temperature T dependence of the 
in-plane resistivity pab for applied magnetic fields H=Q and 
14 T and an applied current 1=1 mA. (b) T dependence 
of magetoresistivity t^pj p where Ap = p(14T) — p(0). (c) T 
dependence of ZFC and PC (H = 50 Oe) magnetization M 
curves, (d) AA//A/ vs. T where AA/ = A/(FC) - Af(ZFC). 
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FIG. 3. (Color online) (a) T dependence of the reduced 
resistivity p{T)/pAip measured at ii" = 0, 0.5, 2, 4, 6, 8, 12 T 
[H II c) and J = 3 mA. {h) H = Q and 6 T for both H \\ c and 
H Jl c with I = 3mA. (c) T and H dependence of the peak 
position for / = 1 and 3 mA. Inset: p,oHc2 — T phase diagram 
of CaFe1.94Coo.06 As2. (From Ref.^^). The solid black lines are 
fits to the GL expression for upper critical field; the dotted 
red lines are fits by WHH relation. 
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ture resistivity normalized by pdip, the value at the lo- 
cal minimum at ^ 13.5 K. Below this temperature the 
resistivity exhibits a local maximum around 10 K, an 
observation that is consistent with other reports. '^'^'■^^ In- 
creasing the magnetic field reveals a semiconducting like 
background with an enhanced resistivity below 10 K. Be- 
low this temperature, however, the resistivity exhibits 
a field-dependent peak and downturn at lower temper- 
ature. The peak position Tpcak (marked with an arrow 
in Figs. 3(a) and 3(b)) shifts to lower T with increasing 
field. Tpcak is suppressed with field and p no longer ex- 
hibits a peak by iJ > 12 T, but continues to rise down to 
the lowest temperature measured. The resistivity peak is 
similar for both H \\ c and H ± c (Fig. 3(b)) indicating 
a very small anisotropy consistent with earlier reports 
in SC samples of this family.^^ Tpeak shifts slightly as a 
function of the applied field direction, and is smaller for 
H\\c. 

The field dependence of Tpcak is summarized in Fig. 
3(c) and bears a striking resemblance to the reported 
upper critical field Hc2 vs T phase diagram for SC 



CaFei.94Coo.o6As2 (upper inset). •^■^ We find H is nearly 
linear in Tpoak with a slope of -2.4 T/K for both cur- 
rent values. Linear extrapolations to zero temperature 
yields critical field values of 18.3 and 16.7 T, respec- 
tively, representing the applied fields necessary to com- 
pletely suppress the downturn in p. In fact, our ob- 
servations are consistent with the development of fila- 
mentary superconductivity, in which only a small vol- 
ume fraction of the material becomes superconducting. 
This fraction is responsible for the suppression of resis- 
tivity, but is not sufficiently extended spatially to lead to 
a diamagnetic signature or a completely zero resistance 
state. The facts that Tpcak is on the same scale as Tc 
in fully superconducting samples and that both temper- 
atures are suppressed with field in similar manners sug- 
gest that the phenomena we observe are not associated 
with an impurity phase. This observation is supported 
by fits to the Ginzburg-Landau (GL) expression for the 
upper critical field, i/c2(T) = i?c2(0)[l - (T/Tc)2]/[1 + 
(T/T^)'^] (solid black lines, yielding ifc2(0) = 15.1 T, 
Tc(0) = 8 K for 1 mA and 12.6 T and 7.6 K for 3 
mA) and the Werthamer-Helfand-Hohenberg (WHH) re- 



lation, Hc2i0) = -0.7Tc{dHc2/dTa) (dotted red lines, 
yielding 13.3 T and 12 T for 1 mA and 3 mA, respec- 
tively). The value of i?c2(0) obtained from GL is larger 
than WHH, a behavior similar to that reported for SC 
CaFei.g4Coo.o6As2. Furthermore, the high temperature 
tail (red dashed lines) in our H — T curve is typical 
for Hc2{T) in the iron-based superconductors.'^^ The fact 
that Tpoak is suppressed with increasing current density 
further suggests that the superconducting filaments are 
weakly pinned. 

It is worth noting that this low temperature anomaly 
around 10 K is about the same for the maximum Tc of 
the material with applied pressure. ^^ Furthermore, it is 
present not only in CaFe2 As2 , but also in the other mem- 
bers of the AFe2As2 family. For example, BaFe2As2 at 
ambient pressure displays a broad maximum in the resis- 
tivity around 20 K, and pressure studies show that even 
very moderate uniaxial stress can induce at least filamen- 
tary SC and a maximum Tc around that temperature.'^^ 
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FIG. 4. (Color online) (a) The resonance frequency of the 
upper satellite of the As along the b direction upon crossing 
an antiphase DW oriented along the a direction, assuming a 
DW thickness S — 5a. (b) The expected lineshape of the As 
resonance for such a DW. 



D. Summary 

Enhanced superfluid density at DWs may be a natural 
consequence of the coupling between SC and AFM 
orders. Indeed, model calculations exhibit an enhanced 
local density of states and superconducting order at 
twin boundaries in the iron pnictides.'^^ It is therefore 
not surprising that superconductivity could nucleate at 
antiphase domain boundaries. However, the fact that 
the emergence of SC coincides with the freezing of the 
AFM DWs is striking, and implies that the former is 
driven by the latter. In other words, defects in the AFM 
background are pinned by the emergence of SC order, 
which would be unstable if the DWs were mobile. 
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by the vector sum of contributions from the four nearest 
neighbor Fe sites. The hyperfine coupling is given by: 



-Hm = jhi • ^ 1. • S(r,), 



(1) 



I. APPENDIX 

A. NMR Spectra of Domain Walls 

In the presence of either a twin boundary DW or an 
antiphase boundary DW, the local magnetic order of the 
Fe moments will be perturbed from the equilibrium anti- 
ferromagnetic structure. This perturbations will modify 
the local hyperfine field at the As sites, which is given 



where the sum is over the four nearest neighbor Fe spins 
S(ri), and the components of the hyperfine tensor B are 
given for CaFe2As2 by: Baa == Bbb = 5.8 kOe//iB, B^c = 
6.0 kOe/nB and Bac = 8.2 kOc///B. 27,37 

The exact structure of the ordered Fe moments in the 
vicinity of the domain walls depends on details of the 
magnetic model and the exchange couplings, which has 
not been completely resolved. However, in general one 
can expect the magnetic order to recover to the equi- 
librium structure within several lattice constants. As a 




FIG. 6. (Color online) Echo decay curves for the As cen- 
tral transition in a field of 5 T for the external field aligned 
perpendicular to the c axis. The echo integral is shown as a 
function of pulse spacing, r, between the 90° and 180° pulses 
at several temperatures, offset vertically for clarity. The tem- 
peratures are (in order of increasing vertical offset): 4.5 K, 6 
K, 8 K, 12.5 K, 15 K, 17.5 K, 20 K, 25 K, and 30 K. The solid 
lines are fits to as described in the text. 



where A is the domain size. This distribution, however, 
should be convoluted with a Gaussian in order to repro- 
duce the experimental data. The width of this Gaussian 
is determined in part by the intrinsic width of the spec- 
trum, as well as the disorder of the domain structure. In 
this case, the spectrum is given by: 
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We fit the experimental spectra to this function with luq, 
a and 6/ ft as variable parameters. We find that cr w 62 
kHz, about 3.5 times the paramagnetic state width. 

The relative weight contained in the low frequency tail 
from the domain walls is given roughly by / = 26/ A. 
For the case of twin domains, transmission electron mi- 
croscopy and optical measurements find domain sizes on 
the order of 40 - 400 nm.^^ Our observation of / « 0.2 
implies that i5 is on the order of 8 - 80 nm, or lOao - 
lOOao, where ao = 0.554 nm is the orthorhombic unit 
cell length. It is important to note that at a twin bound- 
ary, the structural distortion can be fairly small, but the 
magnetic response can be extended spatially depending 
on the details of the magnetic couplings. 



result, the As hyperfine field will have a similar recov- 
ery length at which point the hyperfine field will point 
alternately along ±c. Within the distance 5 from the 
boundary, H^yp will be tilted away from the c-axis and 
have components in the a6-plane. This arrangement 
will modify the resonance frequency, which is propor- 
tional to the vector sum |Ho -I- Hhyp|. We therefore ap- 
proximate the spatial dependence of the resonance fre- 
quency along a direction perpendicular to the DW as 
u!{x) ~ jHq + ujotanh{x/5) where S is the DW width, 
wo ^ 7-ffhyp = 19.2 MHz, and -fHo = 65.4 MHz. For 
example, this would be the case for an antiphase DW 
oriented along the ab direction, where x points along the 
6-direction (see Fig. 1(e)). This dependence is shown in 
Fig. 4(a). For such a one-dimensional model, the his- 
togram of resonance frequencies is shown in Fig. 4(b). 

In order to model the spectrum, we first calculate the 
histogram: 
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5 (lu — ^(a;)) dx, 



(2) 



B. Lineshape in Paramagnetic State 

Fig. 5 shows the spectrum of the upper satellite above 
T/v, which shows a Gaussian form with second moment 
18(1) kHz. 



C. Spin Echo Decay Curves 

The echo decay is shown in Fig. 6 for a scries of 
temperatures, as well as fits to the function A/(2t) = 
Mq exp(— 2r/T2). The echo decay clearly shows exponen- 
tial behavior, and the rate T2 increases with decreasing 
temperature. The temperature dependence is shown in 
Fig. 1(d). 
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